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基于 投影 域 小 波 滤 波 处 理 的 CT 图 像 环 形 伪 影 去 除 方法 
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摘要 :目的 平板 CT 成 像 常 产 生 环形 伪 影 ,其 原因 在 于 平板 探测 器 单元 因 多 种 不 利 因素 的 影响 ,性 能 存在 不 一 致 现象 ,表现 
为 投影 数据 内 含有 直线 结构 信息 。 为 了 实现 平板 CT 优质 成 像 ,本文 提 出 一 种 基于 投影 域 小 波 滤波 处 理 的 环形 伪 影 去 除 方 
法 。 方 法 新 方法 首先 把 线 积分 投影 数据 等 间隔 分 成 四 个 子 投影 数据 ,接着 分 别 对 四 个 子 投影 数据 进行 小 波 分 解 ,生成 小 波 
系数 窍 阵 , 然 后 分 别 对 小 波 系数 进行 滤波 处 理 , 新 的 小 波 系数 矩阵 进行 小 波 重 构 生 成 四 个 新 的 子 投影 数据 ,并 合成 为 新 的 线 
积分 投影 数据 ,最 后 再 用 加 权 均 值 滤波 对 新 的 线 积分 投影 数据 进行 滤波 处 理 ,处 理 后 的 线 积分 投影 数据 采用 滤波 反 投影 

法 进行 图 像 重 建 。 结 果 实验 证 明 测试 算法 有 效 可 行 。 结 论 实验 结果 表明 ,该 方法 能 有 效 地 去 除 环 形 伪 影 ,使 图 像 细 节 清 


晰 可 见 。 
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CT ring artifact reduction using an improved wavelet filtering in the sinogram domain 
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Abstract: Objective Ring artifacts often appear in flat-panel detector-based CT images due to the malfunction or 
mis-calibration of the detector elements that result in stripe artifacts in the line integral projection (sinogram) data. The ring 


artifacts lower the image quality and affect image-based diagnoses. Here we propose a ring artifacts removal approach based 


on wavelet filtering in the sinogram domain. The line integral projection (sinogram) dataset were divided into 4 sub-sinogram 
dataset, and for each of them the wavelet decomposition operation was employed to produce the associated wavelet dataset, 
followed by filtering operation on the vertical detail band and the low-pass detail band. Wavelet reconstruction operation was 
then performed, and the weighted moving average filter was used to yield the filtered sinogram, which was processed using 
filtered back-projection (FBP) for image reconstruction. The results showed that the proposed approach could effectively 


remove the ring artifacts while preserving the structural information of the image. 
Key words: computed tomography; sinogram data; ring artifacts; wavelet filtering; average filter 


环形 盆 影 是 X- 线 计算 机 断层 成 像 常 见 的 一 种 伪 
影 ,其 产生 的 因素 有 很 多 , 主要 包括 探测 需 的 探测 单元 
存在 故障 ,探测 单元 没有 校正 ,探测 器 上 存在 着 灰尘 和 
杂质 等 “2。 环 形 伪 影 在 CT 图 像 中 主要 表现 为 一 系列 
同心 圆 环 或 者 圆 驳 ,其 在 线 积 分 投影 数据 上 表现 为 竖 直 
方向 上 的 直线 ”“。 环 形 伪 影 在 一 定 程度 上 降低 图 像 质 
量 ,给 图 像 噪声 处 理 及 图 像 分 割 等 后 续 处 理 带 来 很 大 困 
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扰 ,严重 影响 影像 诊断 效果 。 因 此 ,去 除 或 者 最 大 程度 
地 减少 环形 伪 影 是 十 分 必要 的 。 

目前 ,大 量 的 环形 伪 影 去 除 方法 相继 被 提出 
主要 分 为 两 类 ,一 类 是 基于 线 积分 投影 数据 的 去 除 方 
法 , 即 图像 重 建 前 的 处 理 方法 。 例 如 ,Boin 等 “采用 均 
值 滤波 方法 去 除 原 始 数 据 中 因 探 测 圳 响应 不 一 致 市 来 
的 高 频 信 号 ,该 方法 可 以 有 效 抑制 环形 伪 影 ,但 图 像 分 
辨 率 会 略 有 下 降 ; 余 晓 钼 等 "采用 Hough 变 换 对 线 积 4 
投影 数据 中 的 环形 伪 影 进行 准确 定位 ,继而 对 相应 位 置 
上 的 数据 进行 线性 插值 实现 校正 ;Miinch 等 “对 线 积 4 
投影 数据 进行 小 波 变换 ,利用 低 通 滤波 方法 滤 除 垂直 方 
向 的 高 频 分 量 , 可 以 实现 弱 强 度 环形 伪 影 的 去 除 。 第 二 
类 是 基于 CT 重建 图 像 的 去 除 方法 ,是 一 种 图 像 后 处 理 
方法 **"。 该 类 方法 的 主要 思想 是 ,将 CT 图 像 通过 坐 
标 变换 ,将 图 像 从 直角 坐标 系 变 换 到 极 坐标 ,在 极 坐标 
下 对 变换 后 的 数据 进行 滤波 处 理 , 去 除 条 形 伪 影 ,最 后 


T 


201712.00798v1 


chinaXiv 


http://www.j-smu.com 


J South Med Univ 2015, 35(9): 1258-1262 


ChinaXiv 合 作 期 刊 
. 1259. 


将 处 理 后 的 数据 再 经 过 坐标 变换 回 到 直角 坐标 系 。 此 
类 方法 需 经 过 两 次 坐标 变换 , 知 滤波 融 的 冰 值 设置 不 习 
当 , 将 导致 图 像 对 比 度 和 分 辨 率 的 明显 下 降 。 

针对 上 述 两 类 方法 各 自 的 不 足 ,本文 提出 一 种 
基于 投影 域 小 波 变换 的 环形 伪 影 去 除 方法 ,此 方法 
把 线 积分 投影 数据 等 间隔 分 成 4 个 子 投影 数据 ,以 
消除 相 邻 探测 带 间 的 影响 ;接着 分 别 对 4 个 子 投影 
数据 进行 小 波 分 解 ,然后 分 别 对 近似 分 量 和 垂直 细 
节 分 量 进行 滤波 处 理 ,以 消除 具有 高 频 特 性 的 直线 
言 息 ;随后 对 处 理 后 的 小 波 系数 矩阵 进行 小 波 重 构 
生成 4 个 新 的 子 投影 数据 ,继而 合成 为 新 的 线 积 4 
投影 数据 ;最 后 再 用 加 权 均 值 滤波 方法 对 新 的 线 积 
分 投影 数据 进行 滤波 处 理 , 人 处 理 后 的 线 积分 投影 数 
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1.1 子 投影 数据 获取 

相 邻 探测 器 通道 的 数据 相关 性 会 影响 条 形 伪 影 
除 效 果 " ,本 文 将 经 对 数 变换 后 的 线 积分 投影 数据 $S 等 
间隔 分 成 上 -个子 投 影 数据 S,, 以 消除 相 邻 探测 需 通 道 的 
数据 相关 性 。 子 投影 数据 获取 过 程 如 下 : 

S (njn)=S(n,7): for 1 <j,, <ja and 

mod(j,L)= mod(i,D) 

其 中 1 <n= Mota < Sjioar, mo 和 juu 分 别 表示 
投影 个 数 和 探测 器 通道 个 数 ,并 且 mod(, 羽 表示; 
(i 和 [做 除法 运算 后 的 余数 。 在 本 文中 ,L=44。 


(1) 


据 采用 滤波 反 投 影 方法 进行 图 像 重建 。 实 验 结 
表明 ,本 文 所 提出 的 方法 可 有 效 地 去 除 环 形 伪 影 。 


1 方法 

基于 小 波 变换 的 CT 图 像 环 形 伪 影 去 除 方法 主要 
包括 3 个 步骤 :(H) 获 取 子 投影 数据 :将 线 积分 投影 数据 
等 间隔 分 成 四 个 子 投影 数据 ;(2) 去 除 条 形 伪 影 :对 每 个 
子 投影 数据 进行 小 波 分 解 ,生成 小 波 系数 矩阵 ,然后 对 
小 波 系数 中 近似 分 量 和 垂直 方向 细节 分 量 进 行 滤波 处 
理 ,新 的 小 波 系数 矩阵 进行 小 波 重 构 生 成 四 个 新 的 子 投 
影 数据 (如 图 1 所 示 );(3) 去 除 残 余 条 形 伪 影 :采用 加 权 
均值 滤波 方法 去 除 残 余 的 条 形 伪 影 以 获得 最 终 线 积分 
投影 数据 。 具 体 步 又 详 述 如 下 。 


图 1 条 形 伪 影 去 除 算 
法 流程 图 

Fig.l Flowchart of the 
procedure to remove 
stripe artifacts. 
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1.2 条 形 伪 影 去 除 

线 积分 投影 数据 经 过 小 波 分 解 后 ,条 形 伪 影 数据 信 
息 主 要 集中 在 近似 分 量 和 垂直 方向 细节 分 量 中 。 
此 ,步骤 2.1 所 获得 的 子 投影 数据 S (i=1,2,3,4) 分 别 进 
行 一 级 小 波 分 解 , 可 以 获得 近似 分 量 ( 如 Aj,i= 1,2,3,4)， 
水 平方 向 细节 分 量 (如 ,i=1,2,3,4) ,垂直 方向 细节 分 
量 ( 如 多 ,=12,3,4) 和 对 角 线 方向 细节 分 量 (如 Di， 
i= 1,2,3,4)( 如 图 1 所 示 )。 对 于 垂直 方向 细节 分 量 的 系 
数 ,先进 行 傅 里 叶 变 换 , 然 后 采用 带 阻尼 滤波 器 处 理 环 
形 伪 影 数据 * 中 。 对 于 近似 分 量 的 系数 ,采用 加 权 均 值 
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滤波 去 除 环形 伪 影 数据 。 计 算 步 又 如 下 : 
(1) 沿 投影 角度 方向 求 和 : 
ai(m) = D>), Ai 站 (2) 
其 中 K,m 表 示 近 似 分 量 维度 大 小 。 

(2) 党 探测 需 单 元 方向 加 权 均 值 滤波 : 


全 加 Ai(m) 

cut A (3) 
GM) = am + -aim), Ovy<0.5 (4) 
ww 下 -总 

te ty (5) 


其 中 6,(m) 表 示 局 部 均值 滤波 结果 ,y 是 一 个 调节 
子 , 在 本 实验 中 ,y 取 值 为 0.3。 

滤波 后 的 小 波 系 数 和 矩阵 进行 小 波 重 构 生 成 新 的 子 
投影 数据 S$;, 进 而 4 个 子 投影 数据 可 以 合成 为 新 的 线 积 
分 投影 数据 $。 
1.3 残余 条 形 伪 影 去 除 

如 果 2.2 中 得 到 的 线 积 分 投影 数据 9 仍然 含有 残余 


图 2 带 有 环形 伪 


影 的 CT 测试 图 像 一 


的 条 形 伪 影 ,我们 可 以 采用 2.2 中 的 加 权 均 值 滤波 方法 
对 线 积分 投影 数据 进行 处 理 以 获得 最 终 处 理 结果 S。 
对 于 最 终 获 取 的 线 积分 投影 数据 $, 采 用 滤波 反 投 影 方 
法 进行 图 像 重 建 。 


2 结果 与 讨论 

为 了 验证 本 文 方法 的 有 效 性 ,我 们 通过 真实 体 模 和 
小 鼠 实 验 分 别 进行 了 定性 分 析 和 定量 分 析 。 实 验 数据 
获取 采用 Siemens Inveon PET-CT 扫 描 设 备 采 集 , 管 电 
压 为 80 kVp, 管 电流 为 100 mAs。 我 们 分 别 对 含有 环 
伪 影 的 真实 体 模 CT 图 像 和 小 鼠 CT 图 像 (图 2) 进 行 环 
上 伪 影 去 除 实验 。 由 于 本 文 算法 与 Munch 提 出 的 方法 
在 环形 伪 影 的 处 理 策略 上 相似 ,我 们 将 对 两 者 进行 试验 
对 比 。 实 验 平台 为 Windows 7, Intel(R) Core (TM) 
i5-2400 处 理 器 ,3.10 GHz,6 GB 内 存 ,程序 设计 采用 
Matlab 编 程 实现 。 


NN 


S| 


Fig.2 Original images with artifacts reconstructed from the raw sinogram 
data. A: Resolution phantom; B: Rat chest. Both of the images are displayed in 


the same window: [0.030 0.055]. 


2.1 体 模 CT 图 像 研究 

图 3A 是 通过 Munch 方 法 校正 后 的 CT 图 像 ,图 3B 
是 通过 本 文 方法 校正 后 的 CT 图 像 。 结 合 Munch 方 法 
特点 , 易 见 ,对 于 强度 较 大 的 环形 伪 影 ,Munch 方 法 难以 
很 好 的 去 除 。 为 了 更 清晰 地 比较 两 种 方法 对 环形 伪 影 
去 除 效果 ,图 4 中 给 出 图 3A 中 所 示 的 感 兴趣 区 域 对 应 
的 局 部 放大 图 。 不 难看 出 ,本文 方 法 能 有 效 的 进行 环形 
伪 影 去 除 , 且 能 保持 图 像 细 节 。 

为 了 进一步 定量 分 析 本 文 方法 ,我 们 计算 了 图 像 的 
局 部 信 噪 比 LLSNR ) ,其 中 计算 公式 如 下 : 

0 
Gan) 
ISNR = = (6) 
0 


0 2 
| 3S P70) 


m=1 m= 


其 中 4 表示 校正 后 图 像 ,0O 表 示 感 兴趣 区 域 (ROTD) 内 
像素 点 个 数 。 其 中 感 兴趣 区 域 如 图 2A 标 记 所 示 。 由 表 
1 的 分 析 可 以 看 出 ,本 文 方法 较 Munch 方 法 提升 了 图 像 
的 局 部 信 噪 比 。 

2.2 小 筷 CT 图 像 研 究 

环形 伪 影 去 除 结果 如 图 5~6 所 示 , 其 中 图 5A 是 小 
鼠 CT 原 始 投 影 数据 ,直线 条 结构 信息 明显 。 图 5B 是 通 
过 Munch 方 法 校正 后 的 CT 投影 数据 ,强度 较 大 的 条 形 
伪 影 依然 存在 。 图 5C 是 通过 本 文 方法 校正 后 的 CT 投 
影 数 据 ,条 形 伪 影 可 以 基本 上 去 除 。 图 6A 是 通过 
Munch 方 法 校正 后 的 CT 图 像 ,图 6B 是 通过 本 文 方法 校 
正 后 的 CT 图 像 。 可 以 看 到 ,我 们 的 方法 可 以 去 除 以 图 
像 中 心 为 圆心 的 环形 伪 影 ,周围 正常 图 像 基本 不 受 影 
响 。 而 Munch 方 法 在 靠近 圆心 的 环形 伪 影 能 够 被 去 
除 , 但 远离 圆心 的 环形 伪 影 无 法 消除 。 
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图 3 体 模 CT 图 像 环形 伪 影 去 除 

Fig.3 Real phantom ring artifact removal study. A: Corrected image by 
Munch method; B: Corrected images by the proposed method. Both of the 
images are displayed in the same window: [0.030 0.055]. 


图 4 对 应 于 图 3 的 不 同方 法 重建 图 像 的 感 兴趣 区 域 
(ROI) 的 局 部 放大 图 


Fig.4 Zoomed ROI images of Fig.3. A: Corrected ee 

image by Munch method; B: Corrected image by the 图 5 小 鼠 CT 投 影 数 据 条 形 伪 影 去 除 

proposed method. Both of the images are displayed Fig.5 Stripe artifact removal from rat chest sinogram 
in the same window: [0.030 0.055]. data. 4: Original sinogram; B: Corrected sinogram by 


Munch method; C: Corrected sinogram by the proposed 
method. Both of the images are displayed in the same 


表 1 不 同方 法 感 兴趣 区 域内 的 局 部 信 品 比 We 


Tab.1 LSNRs of the ROI from different methods 


环形 伪 影 严重 影响 CT 图 像 质量 ,给 图 像 噪声 处 理 


Method Uncorrected Munch method Proposed method 
Oe a 及 图 像 分 割 等 后 续 处 理 带 来 很 大 困扰 ,进而 严重 影响 临 
ROI2 8.53 9.35 9.56 床 医生 对 CT 图 像 的 正确 诊 分 新 。 Miinch 等 5 对 线 积 1 分 投 


影 数据 进 行 小 波 变 换 , 利 用 低 通 滤波 方法 滤 除 垂直 方向 


图 6 小 鼠 CT 图 像 环 形 伪 影 去 除 

Fig.6 Ring artifact removal from rat chest data. 4: Corrected image by Munch 
method; B: Corrected image by the proposed method. Both of the images are 
displayed in the same window: [0.030 0.055]. 
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的 高 频 分 量 ,可 以 实现 弱 强 度 环形 伪 影 的 去 除 ,但 是 对 
于 强度 较 大 的 环形 伪 影 ,Munch 方 法 处 理 的 效果 比较 
差 ,主要 是 因为 在 小 波 变 换 后 ,强度 较 大 的 条 形 伪 影 会 
被 当成 组 织 结构 信息 ,因此 Munch 方 法 不 能 够 消除 这 
些 伪 影 。 本 文 讨 论 分 析 了 Munch 等 提出 的 环形 伪 影 

除 方 法 ,提出 了 一 种 基于 投影 域 小 波 变 换 的 CT 图 像 环 
形 伪 影 去 除 方法 。 新 方法 把 线 积分 投影 数据 等 间隔 分 
成 四 个 子 投影 数据 ,接着 分 别 对 四 个 子 投影 数据 进行 小 
波 分 解 ,然后 分 别 对 近似 分 量 和 垂直 细节 分 量 进 行 滤波 
处 理 , 对 新 的 小 波 系 数 和 矩阵 进行 小 波 重 构 生 成 四 个 新 的 
子 投影 数据 ,并 合成 为 新 的 线 积分 投影 数据 ,最 后 再 用 
加 权 均 值 滤波 对 新 的 线 积 分 投影 数据 进行 滤波 处 理 ,处 
理 后 的 线 积 分 投影 数据 采用 滤波 反 投 影 方法 进行 图 像 
重建 。 基 于 投影 域 小 波 滤波 处 理 的 CT 图 像 环形 伪 影 

除 方法 是 一 种 基于 线 积 分 投影 数据 的 去 除 方法 ,此 方法 
优 于 基于 CT 重建 图 像 的 去 除 方法 **" ,主要 原因 是 经 
过 小 波 滤波 处 理 后 的 投影 数据 通过 滤波 反 投影 方法 进 
行 图 像 重 建 时 ,由 于 投影 数据 中 几乎 不 含有 条 形 伪 影 数 
据 ,使 得 重建 后 的 CT 图 像 的 基本 上 去 除了 环形 伪 影 ,而 
且 能 很 好 地 图 像 分 辨 率 。 试 验 结果 表明 ,Munch 方 法 虽 
然 去 除了 部 分 环形 伪 影 ,但 不 能 去 除 强度 较 大 的 环形 伪 
影 , 而 本 文 提出 的 方法 相 比 Munch 方 法 不 仅 能 很 好 地 
去 除 环形 伪 影 ,还 提高 了 局 部 信 品 比 , 较 好 的 保持 图 像 细 
节 信 息 ,使 图 像 质量 提高 ,有 利于 后 续 处 理 和 量化 分 析 。 
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